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1. Introduction 

Collective modes of motion play a significant role in the low-energy structure of atomic 
nuclei. To account for large quadrupole moments, the spherical shell-model picture 
needed to be extended to spheroidal deformations [1], due to collective polarization 
effects induced by single particles moving in the nuclear medium. As a consequence the 
nucleus can no longer be regarded as a rigid body, rather a soft object with a surface 
that can undergo oscillations and rotations in the laboratory framework. The quantized 
treatment of these excitations in the intrinsic framework [2,3] led to the development 
of the Bohr Hamiltonian in terms of the intrinsic collective variables (3 and 7 [4,5], 
corresponding respectively to the degree of axial and triaxial deformation. The dynamics 
is determined by the potential contained in the Bohr Hamiltonian, which can either 
be constructed from a microscopic theory or through phenomenological considerations. 
When the latter strategy is followed, one can either choose an analytically solvable 
potential, a topic which has recently gained a considerable amount of interest because 
of its application to critical points in phase shape transitions [6-8], or a more general 
expression [4] in terms of the collective variables, determining the surface (CQ). For an 
overview on (approximative) analytically solvable potentials, we would like to refer the 
reader to [9]. 

Analytically solvable potentials are intended as benchmarks in order to study more 
general and complex potentials. To handle these potentials, one needs to perform 
a diagonalization, for which a suitable basis is needed. Within the literature, 
several methods have been proposed and profoundly discussed. Pioneering work 
has been carried out by Bes [10], who determined the explicit 7-soft wavefunctions 
through a coupled differential equation method. Unfortunately, this technique becomes 
tremendously complicated for spin states, higher than L = 6. Therefore, other 
techniques have been developed, fully exploiting the SU(1, 1) x 0(5) structure of the 
five-dimensional harmonic oscillator. Nevertheless, complications arise. Since the 
Hamiltonian is an angular momentum scalar, the eigenstates automatically possess 
good quantum numbers L and M of the physical 0(3) D 0(2) chain which does not 
evolve naturally from the Cartan-Weyl group reduction. As a consequence one is forced 
to construct explicit wavefunctions, starting either from basic building blocks [11-13] 
or from a projective coherent state formalism [14], constituting an orthonormal basis 
[15,16]. Nevertheless, the alternative Cartan-Weyl [17,18] reduction path may also be 
followed as it leads to a reduction scheme which is more natural in a mathematical 
sense, though the physical meaning of the quantum numbers is partially lost. This 
strategy was followed by Hecht [19] to construct fractional parentage coefficients for 
spin-2 phonons, that were used by the Frankfurt group [4, 20] in the development of the 
General Collective model. 

It is noteworthy that new techniques have been proposed within the last decennium. 
First, the vector coherent state formalism [21-24] and much more recently the algebraic 
tractable model [25, 26] were developed, enabling the construction of the quadrupole 
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harmonic oscillator representations, exhibiting good 0(3) angular momentum quantum 
numbers. 

In the present paper, the path of the natural Cartan-Weyl reduction is taken. It will 
be shown that the matrix elements of the quadrupole variable can be extracted within 
this basis, without making use of the explicit representations in terms of the collective 
variables. In the end, it will turn out that the basic commutation relations of the 
collective variables suffice to fix the complete structure of the algebra, and furthermore 
the dynamics of the Hamiltonian. 



2. The collective model 



Within the framework of the geometrical model, the nucleus is regarded as a liquid drop 
with a surface R(9, 0) described by a multipole expansion using spherical harmonics in 
the laboratory system 

R(6, ^ = R [l + J2 <t>)) > (!) 

A, fi 

which defines the set of collective coordinates a\ tfl of multipolarity A and projection 
\i. Up to quadrupole deformation, the surface JTJ is restricted to spheroidal deforma- 
tions determined by the variables a2, M which will be abbreviated to a M from here on. 
Although the intrinsic surface is unambiguously described by this set of variables, it is 
convenient to rotate from the laboratory to the intrinsic framework by means of the 
Euler angles (9i,i = 1,2,3). Doing so, the collective variables (/?, 7) are introduced as 
intrinsic parameters of the ellipsoid, rendering a straightforward interpretation of axial 
and triaxial deformation [4]. 



a. 



(2) 



This set of collective variables is sufficient for the determination of the static properties 
of nuclear shapes. To build in the essential quantum mechanical dynamics, canonic 
conjugate momenta 7iy need to be incorporated. These must fulfill the standard 
commutation relations [4] 

[V> a d = -^<W> [V' 7 ^]^ ) [cy,a M ]=0. (3) 

Note that the variables have become operators though we silently omit the operator 
symbol to avoid notational overload. 

To establish the 577(1,1) x 0(5) group structure, it is convenient to introduce the 
following recoupling formula 

(a ■ a)(ir* ■ ir*) =(a ■ ir*)(a ■ tt*) + 3ih(a ■ 71*) 



- 2([an 



*l(i) 



*l(i) 



+ [air' 



1(3) 



*i(3h 



(4) 



where the complex conjugate ir* is introduced to ensure for good angular momentum 
transformation properties and 2] ■ Si — (— ) / v / 2T+T[T/S , i]o°' ) 
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The 3 operators (a ■ a), (n* ■ rr*) and (a ■ 7r*) generate the algebra of an SU(1, 1) group, 
which forms a direct product together with the 0(5) group, built from the 10 operators 
[ck7T*] ^ and [tt7T*]^. Whereas the SU (1,1) group is strongly linked to the excitations 
in the radial variable f3, the 0(5) group encompasses the 7 vibrations coupled to the 
rotational structure. In this work, we concentrate on the application of the Cartan-Weyl 
scheme on the 0(5) group, leaving a freedom of choice of a suitable SU(1, 1) basis. 



3. The Cartan-Weyl reduction of 0(5) 

The commutation relations of the operators L M and Om, defined by 

r„ //7r *i (1) _ Jh_ j 



= 7M° 



(5) 



'Mi 



span the algebra of the 0(5) group. 

[L m , L m >] = -V2(lmlm'\lm + m)L m+m <, (6) 

[L m , O m i] = -2\/3(lm3m'|3m + m')O m+m >, (7) 
[O m , O m >] = —2V7 (3m3m'\lm + m')L m+m < 

+ \/6{3m3m'\3m + m')O m+m i. (8) 

The M = projections {L , O } form an intuitive choice of the Cartan subalgebra within 
the set {L m , O m >}. This set has the advantage of incorporating the angular momentum 
projection operator L in the physical group reduction chain 0(5) D 0(3) D 0(2). 
Nevertheless, it is not explicitly contained in the natural Cartan-Weyl reduction chain 
0(5) D 0(4) = SU(2) x SU(2) which can be realized through the following rotation 
[11,21] 

X + = -l(V2L +1 + VSO +1 ), Y + = -^0 +3 , 
X_ = l(V2L_, + >/3CLi), Y_ = ^0_ 3 , 

X = ^(L + 3O ), F = ^(3L -Oo), (9) 

T hh = ^0 +2 , Tii = -4b(V3L +1 - v/20 +1 ), 



22 V10 _rz '' -22 V50 

T -\-\ = -TTo ^ T \~\ = Wo^ L -i ~ ^O-i). 

The group reduction is immediately clear, as the sets {X , X±} and {Y , Y±} both span 
standard SU (2) algebras. Furthermore all generators of the one SU(2) algebra commute 
with all generators of the other. The commutation relations are given by 

[X ,X ± ]=±X ± ,[X+,X_] = 2X , 

[Y ,Y ± ] = ±Y ± , [Y + ,Y_]=2Y , (10) 
[X , Y ] = 0, [X ± , Y ± ] = [X ± , Y T ] = 0. 

So the reduction is 0(5) D 0(4) = SU(2) X x 5C/(2) y . The non-0(4) operators T M!y 
can be identified as the 4 components of a bitensor of character {|, |} within the 
SU(2) x S , f/(2) scheme, according to Racah's definition [27]. The index \i denotes 
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the bitensor component relative to the SU(2)x group, while v is the component with 
respect to SU(2) Y 

[Xq, T^v] = fiT^, 

[X ± ,T M „] = ^(i-F^di/x+l)^, ^ 

[Y ± ,T, U ] = j{\^v){\±v + l)T ta/±1 . 

The internal commutation relations of the T bitensor completes the Cartan-Weyl 
structure, which can be found in table (pQ) 



* 


T i i 

9. 2 


Ti i 

2 2 


Tii Tn 

2 2 2 2 


T i i 

2 2 







gX_ ^(^o + yo) 


Ti i 

2 2 







|(Xo — yo) — 


T ii 

2 2 


-§x_ 


— 2"(^o — yo) 


o -iy + 


Tn 

2 2 


-i(X + Yo) 


|x + 


±y+ o 



Table 1. Multiplication table for the internal commutation relations of the T bitensor. 
The multiplication * symbolizes the standard commutation. 



Once the commutation relations have been determined within the Cartan-Weyl basis, 
it is instructive to construct the root diagram. Figure Q] shows 2 different realizations of 
the same root diagram, depending on the choice of the Cartan subalgebra. On the left 
side (Fig. UK) a standard root diagram with respect to the {X , Y } Cartan subalgebra is 
depicted, while on the right side (Fig. [lb), a more physical subalgebra {L , Oq) is chosen 
as a reference frame. The latter framework has a visual advantage, since the projection 
of the generators on the Lo-axis is readily established. This enhances the insight in the 
problem of constructing wavefunctions with good angular momentum from the weight 
diagrams in the Cartan-Weyl basis (see section (El)). 



4. Representations of 0(5) 

Every subgroup in the group reduction chain provides an associated Casimir operator. 
The quadratic Casimir operator of 0(5) can be constructed from the Killing form [18] 

C 2 [0(5)] = |(T-T + 0-0), (12) 
= 2(X 2 + y 2 -2[TT](°°)). (13) 

with [TT](°°) denoting the scalar Clebsch Gordan coupling with respect to both SU{2) X 
and SU(2) Y , and X 2 and Y 2 the quadratic Casimir operator of the respective SU(2) 
groups 

X 2 = X 2 + |(X + X_ + X_X + ), (14) 

y 2 = y 2 + |(y + y- + y_y+). (15) 
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Figure 1. The root diagrams of the 0(5) algebra in the Cartan-Weyl basis for cither 
the (a) natural {X ,lo} or (b) physical {L ,Oq} Cartan subalgebra 



Starting from the explicit expressions of the generators (see Appendix A) in terms of the 



collective variables and the canonic conjugate momenta, the following operator identity 
can be proven 

X 2 -Y 2 = } (16) 

which is true in general for symmetric representations [11]. The consequence of this 
identity is that we are left with 4 operators that commute among each others, i.e. 
the quadratic Casimir operator of 0(5), the quadratic Casimir operator of SU{2) X 
and SU{2) Y (X 2 = Y 2 ) and the Cartan subalgebra {X , Y } which are the respective 
linear Casimir operators of the 0{2) x and 0(2)y subgroups. As a result, we obtain a 
representation which is determined by 4 independent quantum numbers 

\vXM x M Y ) (17) 

with 

C 2 [0(5)}\vXM x M Y ) = v(v + 3)\vXM x M Y ), (18) 

X 2 \vXM x M Y ) = Y 2 \vXM x M Y ) = X(X + l)\vXM x M Y ), (19) 

X \vXM x M Y ) = M x \vXM x M Y ), (20) 

Y \vXM x M Y ) = M Y \vXM x M Y ). (21) 

Now that the basis to work in is fixed, we can study the action of the generators as 
they hop through the representations with fixed quantum number v. Acting with the 
0(4) = SU {2) x x SU (2)y generators on \v XM X M Y ) is trivial because of the well-known 
angular momentum theory 

X±\vXM x M Y ) = y/(X T M X )(X ±M X + l)\vXM x ± 1, My), (22) 
X \vXM x M Y ) = M x \vXM x M Y } } (23) 
Y±\vXM x M Y ) = y/(X =f M Y )(X ± M Y + l)\vXM x ,M Y ± 1), (24) 
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Y \vXM x M Y ) = M Y \vXM x M Y ). (25) 

The action of T^ u on \vX(Mx, My)) is less trivial, though the bitensorial character of 
T can be well exploited. Since T is a {\\} bitensor, it can only connect representations 
that differ \ in quantum number X 

X + i. M v + //, Mv 4- v\ 

(26) 



T^\vXM x M Y ) = a\v,X + §, M x + My + v) 



+ 6|v, X - §, M x + /i, My + u) 

The coefficients a and b are not only dependent on v and X, but also on the projection 
quantum numbers /i, z/, Mx and My. However, these projections can be filtered out 
by means of the Wigner-Eckart theorem. As the SU(2)x forms a direct product with 
SU(2)y, we can apply the theorem for both groups, independently from each other. 
As a result, the dependency on the projection quantum numbers is completely factored 
out in the Wigner-3j symbols. This leaves a double reduced matrix element [j] to be 
calculated. 



a = (v, X + |, M x + fi, M Y + v\T^\vXM x M Y ) 




- V V 

(v, X - |, M x + (i, M Y + u\T^\vXM x M Y ) 




1 (27) 
(vX+-\\T\\vX), 




1 (28) 
with k — 2X + 1 — Mx — My — \i — v. 

In order to calculate the double reduced matrix element, we have 2 types of expressions 
at hand. On the one hand the internal commutation relations of the T bitensor (see 
Tabled]) and on the other hand the Casimir operator of 0(5) (fT3j) . First we consider the 
internal commutation relations, in which case it is instructive to proceed by means of 
an example although the obtained result is generally valid. Take e.g. the commutation 
relation [T_i_i,Tii] = |(Xq + 1q), an d sandwich it with the state \vXMxM Y ) 



2 2 2 2 J 



(vXMxMy\T_i_iTii\vXM x My) - (vXMxMy\TiiT_i_i\vXM x My) 

s 2222 2222 

= \{M X + M Y ). 

At this point, we can insert a complete set of intermediate states \v' X' M' X M Y ) between 
the two generators. 

( vXM x M Y I T_ i _ 1 1 v' X 1 M' x M' Y ) (v'X' M' x M' Y \ Tx 1 1 vXM x M Y ) 

v'X'M' x M' Y 

- (vXMxMy\Tii\v'X'M x My)(v'X'M'xMy\T_i_i\vXMxMy) (30) 

--\{M X + M Y ). 



% We formally use the single reduced matrix notation in order to express the double reduced matrix, 
as any confusion between normal and double reduced matrix element is excluded within this work. 
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Due to symmetry considerations, a large amount of the matrix elements in the 
summation are identically zero. First of all the SU(2) X x SU(2)y bitensor character 
of T dictates strict selection rules with respect to X, M x and My. As a result the 
summation over X', M' x and M' x is restricted to specific values which are completely 
governed by the Wigner-3j symbol in (127|28p . Secondly, the components T^ u of T are 
0(5) generators, which cannot alter the seniority quantum number v. So, the summation 
over v' is reduced to one state v' — v. 

Once the restriction in the summation is carried out, it is convenient to apply the 
Wigner-Eckart theorem (I27II28P and after some tedious algebra we obtain a relationship 
for the double reduced matrix elements 

(vX\\T\\vX+ \){vX + \\\T\\vX) {vX\\T\\vX -\){vX -\\\T\\vX) (2X + l) 2 



.(31) 



2X + 2 2X 2 

The same procedure can be followed for the quadratic Casimir of 0(5). Sandwiching 
equation ([TBI with \vXM x M Y ) yields 

(vXM X My\(T^iiTi_i +Tx_lT_ii -T_i_iTii-TiiT_x_i)\vXM X My) 

s 2 2 2 2 22 22 2222 22 22 (32) 

= v(v + 3) -AX(X + 1). 
By inserting again a complete set, applying the Wigner-Eckart theorem and making use 
of the previously derived relation (13"T1) . we obtain the result 

-4(vX\\T\\vX+±)(vX + l\\T\\vX) = (v - 2X)(v + 2X + 3){2X + 1)(2X + 2). (33) 

This can slightly be rewritten, if one takes the Hermitian conjugate of the T bitensor 
into consideration. 

Tl = {-iy +v T_,_ u . (34) 

It can be proven that this leads towards the following expression for the double reduced 
matrix elements 

{v,X\\T\\v,X+\)* = -(v,X + l\\T\\v,X). (35) 

As a result, we can write 

|(*;X||T||t;X+i>| 2 = \{v - 2X)(v + 2X + 3)(2X + 1)(2X + 2), (36) 

\(vX\\T\\vX - \)\ 2 = \{v -2X + l)(v + 2X + 2)(2X)(2X + 1). (37) 

So the double reduced matrix elements are determined up to a phase. Here we fix the 
relative sign of (vX\ |T| + and (v X\ \T\ \vX--z) to be opposite, as it is the only way 
to obtain eigenstates with real angular momentum L in the physical basis (see section 

©). 

Once that the double reduced matrix elements are determined, they can be plugged into 
equations (I27ll28p . yielding the action of the bitensor T components. 

Tii\vXM x M Y ) = 

"^ X + Mx + l)[X ^My + l)(v . 2X){v+2X + 3) 

2v/(2X + l)(2X + 2) ' 2 2 2/ (38) 

y/(X - M X )(X - M Y )(v -2X + l)(v + 2X + 2) 



2v/(2X)(2X + l) 



\vX-lM x + ±M Y + ±), 
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Ti_i\vXM x M Y ) = 

• = V (X + M X l 1)(X M„ + !)(„ - 2X)(„ + 2X + 3) _ 

2v/(2X + l)(2X + 2) f39) 
■ V(^-Mx)(X + M 1 .)(,-2A- + l)(,. + 2X + 2) |;) _ Y _ ^ + ^ _ 



2v/(2X)(2X + l 
T_ii|t;XM x My) 



- Mx + Igra + ^ - 2X)Q + 2X + 3) 

2^(2X + l)(2X + 2) 2 2 f40) 

V(X + M x )(X-M y) (»,-2X + l)(», + 2X + 2) |t;X _ ^_ y + 



2v/(2X)(2X + l 
T_i_i|wXM x M y ) 



V(X-M X + 1)(X-M Y + 1)( V -2X)( V + 2XT3) 

2v/(2X + l)(2X + 2) f4l) 
y/(X + M X )(X + Af y )(t; - 2X + l)(t; + 2X + 2) x 

2v/(2X)(2X+ 1) ' 2 2 2/ 

From these expressions it is clearly seen that no representations can be constructed 
with X > |, as the representations must have a positive definite norm. Combining 
these results with the standard quantum reduction rules for the SU(2) group, we can 
label all basis states of a representation with fixed v as follows 

X = 0...v/2, 

M x = -X . . . X, (42) 
My = -X ... X 

Figure (T5]) gives a visual interpretation of the reduction rules for the representation 
v = 2. 



5. Matrix elements of collective variables 

The Hamiltonian describing a system undergoing quadrupole collective excitations 
contains a potential V(a), written in terms of the collective variables a^. Indeed, it 
turns out that 

[aa] (2) -a ~/5 3 cos37, (43) 

can be considered as the building block of the 7 part of the potential V(fi, 7) in the 
intrinsic frame [4]. Therefore, matrix elements of within a suitable basis are needed 
for the construction of the matrix representation of the Hamiltonian. As the chosen 
framework in the present paper is the Cartan-Weyl natural basis, we proceed within 
this basis and show that all matrix elements can be calculated by means of an algebraic 
procedure, similar to the one proposed in the preceding section. 

First, we need to establish the bitensor character of the collective variables a M with 
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Figure 2. Visual interpretation of an 0(5) representation with v — 2. Every sphere 
denotes a single basis state. The representation is organized in planes with distinct X 
quantum number, which contain (2X + l) 2 (Mx,My) projection states. 



respect to SU(2)x x SU(2) Y . Calculating the commutation relations of a M with the 
SU (2) generators (which is done most conveniently using the explicit expressions given 
Appendix A| ), we can summerize them as 



m 



[X Q ,a$\ = iw$, (44) 



[X ± , <#] = ^(\T»)(X±» + l)a£ lv , (45) 
[Y ,<] = ^<, (46) 
[Y ± , c$) = v/(A T ,)(A±, + l)aJ ±1 , (47) 

where the 5 collective variables have been relabelled as follows 

ii ii ii ii 

ct 2 = all, ^ = a"n,a-i = aPi,a_ 2 = \, (48) 

2 2 22 22 2 2 J 

00 * 



a = aZ)- (49) 

This clearly states that the 5 projections of a can be divided into the 4 components 
of a {2I} bispinor and a single biscalar, according to Racah [27]. We can again define 
double reduced matrix elements 

(vXM x M Y \a^\v'X'M x M Y ) 

X X ' \ I y\\ All >y>\ (50) 
A/fl (vX\\a \\vX), 

V My J 

with k = 2X — Mx — My- It is noteworthy that, contrary to the matrix elements of the 
generators T^ u , v' is not necessarily equal to v. To obtain explicit expressions for the 
double reduced matrix elements, we start from the commutation relations 




'■1 

j 1 1 —\ 



i T ^i a ¥u>] = — Ty|— 5-w'5-w'a O0 , (51) 
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It 1 



001 
00 J 



-.a 



1 1 

2 2 



r„ AA ^A'A'i 

[<V> °wJ 



= 0. 
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(52) 
(53) 



First a relationship between (t>X||a°||u 'X') and (uXHaaJlu' X') needs to be established. 
This can be accomplished using the commutation relations (I5TT) . We consider the specific 



ii 

case pii,aifi] = and construct the following matrix elements 

22 2 2 

(vX ±\,M X + |, M y + ||[Tii, a\\]\v'X =F |, M x — \, M Y -\), (54) 
(i;XM x My|[Tii,af!]|t/XM x - l,M y - 1), (55) 

2 2 2 2 

in terms of the double reduced matrix elements. This can be achieved by inserting 

1 1 

a complete set of basis states between the generator Tii and the variable off, then 

2 2 2 2 

making use of the T matrix elements (1381) obtained in the previous section and the 
Wigner-Eckart theorem fl50l) . The outcome of these tedious although straightforward 
calculations are, respectively for fl54|) (with minus and positive sign) and fl55l) 

y /(v'-2X)(v' + 2X + 3) 



(vX - l\\a»\\v'X)- 



*2X 



■{vX\\ai\\v l X + \) 



1 ^{v -2X + l){v + 2X + 2) 



(56) 



V2X + 2 



2 V2X + 2 



-(vX\\a*\\v'X-l) 
(vX + ±\\a*\\v'X) 



1 y/(v-2X)(v + 2X + 3) 



(57) 



+ 



2X)(v + 2X + 3) 
(vX\\a^\\v'X + \) 



V2X 

(v' + l)(v' + 2) _ (v + l)(v + 2) 
2X 2X + 2 

(v + l)(v + 2) (t/ + l)(t/ + 2) 



2X 



2X + 2 



y /(v'-2X)(v' + 2X + 3) 

1 1 

The same procedure can be repeated for the commutation relation \Ti i, a 2 \ i 

22 -jj 

obtain again (156"]) and (J57j) . accompanied by the following expression 
(vJf + |||a3||v / X> 



(58) 
= 0. 

0. We 



y/{v-2X){v + 2X + 3) 
(vX\\a^\\v'X+ \) 



[{v - 2X){v + 2X + 3) - (v' -2X + l)(v' + 2X + 2)] 

[{v -2X + l)(v + 2X + 2)- {v' - 2X){v' + 2X + 3)] = 



vV-2X)(w' + 2X + 3) 

Combining (!58|) with (I59I) gives a homogeneous set of two equations in two variables 
(vX + || \az | and (vX\ \a% \ \v'X + |), rendering the trivial zero solution, unless the 
determinant of the associated matrix identically vanishes. This is only possible when 
v' — v ± 1, which proves the common knowledge that a forms an 0(5)-tensor of rank 1. 



]uadrupole collective variables in the natural Cartan-Weyl basis 12 



1 1 



Finally, we repeat the procedure for [Ti i,o 2 i i] — — 75 ^oo' Besides ( 15 6 p and (1571) . we 

2 2 — 5~ 2 » 

obtain the expression 
(vX + |||at||?/X) 



y>-2X)(?; + 2X + 3) 



[X(v' + l)(v' + 2) - (X + l)(u + \){v + 2) + 2(2X + l) 2 ] 



+ ^IMII^ + S) + 1)( , + 2 ) _ (X + !)(,/ + + 2) + 2(2X + l) 2 ] (6 ° } 

y/iv 1 -2X){v' + 2X + 3) L J 



= -2v / 2 v / ( 2X + !)( 2X + 2)(t;X||a ||t/X). 
Solving the set of equations (1311 and (16"Uj) (or equivalently (j5*Sl) and (J5D|) ) results in 
expressions of all possible double reduced matrix elements of a 2 as a function of the 
double reduced matrix elements of a . 

From this point on, we will explicitly take into account that the a variable connects 
representations with Av = 1, omitting the other matrix elements which are identically 
zero. For v' — v + 1 we obtain 



{v ,X + i \\ a inv + l,X) = -y§^^^^ { vXna^ + l,X), (61) 

(v,X\\ai\\v + l,X + I) = -Ly||±iyi±g±l( X||a°|| B + 1,.Y>, (62) 
and for v' = v — 1 



+ ||M||«- = ^V2^V "t-1x 3 < " X||a011 " " ^ <63) 

- 1,A- + |> = --Lyf±f ^HfH| ( .„A-|| a »||. - 1, A'). (64) 

So, we only have to determine the biscalar double reduced matrix elements. Although 
the commutation relations (1531) seem trivial, they are convenient in the derivation of the 
{00} double reduced matrix elements. If we consider only the non-trivial commutation 
relations for which A = | and A' = |, and apply again the same procedure which has 
been used throughout the present paper, we obtain the following result 

E„, (vX\\ah | \ V 'X + \) (v'X + \vX) 

2X + 2 /g^N 

EA^W^Wv'x - \)(v' x - \\\€&\\vx) _ Q 



2X 

Now, taking all derived expressions (I56l) . (l57|) and (ISTj) to (JM]) into account, we can 
rewrite the relation (l65l) as 

(„ - 2X | 2 + 2X + 3) W 1- + + , , 

(2v + l) (66) 
(tX:||a ||?;- IX) (u- 1X| |t;X>. 



(v + 2X + 2)(?;-2X) 

This relation differs from the previously derived expressions with respect to the quantum 
numbers. The expressions ( I56p to (1641) relate matrix elements with different X 
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connections, though the seniority connection (v to v' = v ± 1) was fixed. Now, 
relates matrix elements with different seniority connection, leaving the X quantum 
number unaltered. 

At last, in order to obtain explicit expressions, we return to the geometry of the problem. 
It has been mentioned earlier that the operator a ■ a commutes with all the generators 
of 0(5), making it an 0(5) scalar. Therefore, this operator can be treated as a constant 
with respect to the 0(5) scheme. We call this constant (3 2 , referring to the radial 
deformation parameter in (j2J). As a consequence, we can write 

(vXM x M Y \a ■ a\vXM x M Y ) = f3 2 . (67) 

The procedure used culminates into closed expressions of the matrix elements. By 
inserting a complete set of basis states between the variables of (1671) . we can rewrite this 
expression in terms of double reduced matrix elements. 

1 = (2X X +1)2 5>X| |a°| \v'X) {v'X\\a?\ \vX) 

v' 

+ (vX\\a* | \v'X + \) (v'X + || |a* | \vX) (68) 

+ (vX\\a^\\v'X - \){v'X - \\\afc\\vX)\. 

All these different matrix elements can be reduced to a single one by means of the 
reduction rules (156]) to fl66l) . As a result, we obtain 



(v,X\\oP\\v + l,X)(v + l,X\\a°\\v,X) = ~ + 1)(, + 2X + 3) 

\ , II II > /\ ,1111,/ (2w + 3)(2u + 5) V ; 



(69) 



(v,X\\a \\v-l,X)(v-l,X\\a°\\v,X)= {v 2X J[ V ± 2 * + 2) (2X + I) 2 /? 2 . (70) 

[2v + lj{2v + 3) 

Taking into account that a\ = ocq, we can write 

(vX\\a°\\v'X)* = (v'X\\a°\\vX), (71) 



and summerize 



2X + l)(u + 2X + 3), ov 



| (v, X\ \a"\ \v + 1, X)\ = v - - I ' ; (2A + l) 2 /3 2 , (72) 

(2t> + 3)(2t> + 5) 

1 1 nil Vl9 (v -2X)(v + 2X + 2) , x9 „ 



which is equivalent to 



<„,«„- !,*> = / ^^±M±j) (2 X + 1),, (75) 

y (2t> + l){2v + 3) 

as «o is a hermitian operator. 

Making use of the reduction rules (I61ti64p . we can construct all double reduced matrix 
elements of the a variable. Taking the appropriate Wigner-3j coefficients into account, 
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the total matrix elements of the a variable can easily be derived. As an example we 
evaluate the matrix element 



{vXM x My\o^WXM x My) 




(vX\\a°\\v',X), 



(76) 



(with k = 2X — M x — My). This leads to the closed expression 



(vXM x M Y \a™\v + 1, XM X M Y ) = ^ ' g + ^+5) + ~ ' (?7) 




(vXM x My\c®\v - 1,XM X M Y ) =P,I {V (2 ^ } ^ 2 + 3 | 2) - (78) 

There is a subtlety involved with equation (167j) . /3 2 can be either be regarded as the 
radial variable in the 5-dimensional Euclidean space, which is a constant by definition 
under rotations of the 0(5) orthogonal group, or it can be recognized as a generator 
of the aforementioned SU(1, 1) algebra. In the latter scheme, the 0(5) Hilbert space 
needs to be extended to incorporate this SU(1, 1) basis. Then, it is more convenient to 
move over to a boson creation and annihilation realization 

with [6 M , bl] = and 6 M = (— 1)^6_ M , as it gives immediately rise to the SU(1, 1) 
algebra spanned by 

B + = \b ] -b\ B_ = \b-~b, S = i(6t.6 + 6-6 t ), (80) 

which is closely connected to the Hamiltonian of a spherical harmonic oscillator via 
the Cartan operator Bq. A recent study [26] has proven that an extension to a David- 
son based SU(1, 1) algebra has a major advantage over the spherical basis in the sense 
that the (3 deformation of atomic nuclei is naturally included in the basis representa- 
tions, leading to a higher convergency speed when dealing with actual deformations in 
collective nuclear model calculations. Therefore, we prefer to regard f3 2 as a constant 
over 0(5) rotations, as no bosonic representations of the Davidson SU(1, 1) are known. 
Moreover recent developments [28] in the framework of the factorization method [29] 
have shown that it is possible to extract the needed (3 matrix elements by means of 
generalized raising and lowering operators for a Davidson type of potentials. 



6. Rotation to the physical basis 

The drawback of the Cartan-Weyl reduction is that its basis is not naturally 
compatible with the physical angular momentum quantum number L which emerges 
from experimental energy spectra. This comes from the fact that L-L does not commute 
with X 2 , implying that a basis diagonalizing both operators is non-existent. Therefore, 
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a rotation from the natural group chain 0(5) D 0(4) = SU(2) x 377(2) to the physical 
chain 0(5) D 0(3) D 0(2) is needed. Fortunately, the Casimir operator L associated 
with the physical 0(2) group is diagonal in the Cartan-Weyl basis, leaving L ■ L the 
only operator to diagonalize. 

L-L = L 2 + \(L + L_ + L_L + ). (81) 

Rewriting the 0(3) generators in terms of the Cartan-Weyl generators gives 

L ± = 2X ± + v / 12T T i ± i, (82) 

L = X + 3Y , (83) 

so that L ■ L can be written as 

L-L = 4X 2 -3[(X -3Y + ±)(X + Y + l)-\] 

84 

+ 4 v / 3[T_iiX_ +Ti_iX+l + 12Ti_iT_ii. 

2 2 22 22 22 

The action of all generators are known in the natural basis (see section H]), so the matrix 
elements of a matrix representation can easily be calculated. The dimension of the 
matrix is governed by the seniority quantum number since the generators involved in 
the expression L ■ L cannot alter the 0(5) quantum number, which means that there 
is an associated matrix with every v. This matrix is even further reducible if one 
takes the Lq operator into account. Since Lq can be written as Lq = Xq + 3Yq, it 
is immediately diagonal in the Cartan-Weyl basis, making M = Mx + 3My a good 
quantum number. As a consequence, the total matrix representation of L ■ L can be 
divided in separate sub-matrices with distinct M quantum number. The possible basis 
states \vXMxMy) spanning the sub-matrices with M = Mx + 3My can easily be 
recognized in the tilted weight diagrams (Figure [3]). The diagram is tilted with respect 
to the angular momentum operator L , so that the vertical projection of every basis 
state immediately gives the L component. As a result, all basis states, lying on the 
same vertical projection line form a subspace of states for which M = Mx + 3My holds. 
Although the rotation from the natural towards the physical basis corresponds to a 
standard diagonalization problem, it is useful to study some specific cases. It is readily 
seen from figure [3] that the projection M — 2v can only be constructed from one single 
basis state \v, X — |, M x = |, M Y = f ), as there is only one projection state. Therefore, 
the matrix representation is one dimensional 

(v, !(ff)|£ ■ L\v, f (f f )) = 2v{2v + 1). (85) 

The same is valid for M = 2v — 1. The only basis state with this M projection 
is \v,X = |,Mx = | — I, M Y = |), giving the same eigenvalue 2v(2v + 1). For 
M = 2v - 2, there are two different states: \v, X = |, M x = f - 2, My = f) and 
\v,X = | — \,M X = \ — |, My = I — |). which gives a 2 dimensional matrix with 
eigenvalues and accompanying eigenvectors 



A + = 2„(2„ + l) -(VfeWl^i), (86) 



A_ = (2v - 2){2v - 1) - ( y/^ t -yj&=g) • (87) 
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Figure 3. Projections of the v = 4 irreps on the Lq axis, connecting all basis states 
for which M = Mx + 3My is equal. The line most on the right only projects the 
Mx = v/2, My = v/2 state on the Lq axis, resulting in an M = 2v state. For v = 4 
the maximal Lq projection corresponds to M = 8 

So it is clear that the associated eigenvector of A + belongs to the L — 2v multiplet while 
the eigenvector of A_ will be the heighest M state of the L = 2v — 2 multiplet. Basically, 
this procedure can be repeated up to M = by means of a symbolic mathematical 
computer program or by means of numerical procedures. 

Finally we discuss the dimension of the matrix representations, as they are important 
in actual calculations. The total number of basis states within a representation v can 
be determined as 

v/2 

£(2X + l) 2 = -(2v + 3)(v + 2)(v + 1). (88) 
x=o 

However, the dimensions of the M = sub-blocks are lower than the total v 
representation space. These dimensions can be calculated according to the following 
formula, depending whether v is even or odd. If we define X m = v/2, the number of 
M = projections is then given by 

3(X m |3) 2 + [X m mod3 + l][2(X m |3) + 1] + 3[(X m + 1)|3 - \][{X m + 1)|3] 

(89) 

+ 2[(X m + l)\3][(X m + l)mod3 + 1], 
for even v. For odd v we obtain 

3[(X m + §)|3 - l][(X m + |)|3] + 2[(X m + |)|3][(X m + |)mod3 + 1] 

+ 3[(X m - i)|3] 2 + [(X m - i)mod3 + l][2((X m - |)|3) + 1]. 
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These two dimension formulas are plotted in figure HI From this figure, it is clear that the 



le+06 p 




1 10 100 

V 

Figure 4. Total dimension of a representations with a given seniority v, compared to 
the dimension of the subspace of states with M — 0. It is clear that the dimension 
of the M = subspace increases quadratically with increasing seniority, however it 
remains feasible for realistic calculations. For v = 100, the dimension of M — space 
is 1717, while the total dimension is 348551. 

dimension of the M = subspace stays reasonable with respect to modern computation 
standards, as long as relatively low-order seniorities are considered. Anyhow, when 
performing realistic calculations, the transformation from natural to physical basis 
does not need to be repeated for every calculation, as the rotation is independent of 
the specific physical system (Hamiltonian) under study. In practical calculations, the 
rotation only has to be carried out once and stored for later use. 

7. Conclusions and outlook 

Collective modes of motions have proven to be very important in atomic nuclei, 
away from the shell-closures. Therefore it is of major interest to construct schematic 
Hamiltonians in the significant degrees of freedom. Throughout the last decades much 
attention has been given to solving the Schrodinger equation for general potential energy 
surfaces in the collective variables. This resulted in a number of techniques, based on 
combinations of analytical and algebraic considerations. The present manuscript adds 
a method which is completely algebraic in the sense that no normalized highest weight 
states need to be constructed. As a matter of fact, although the 7- rotational structure 
of the collective model is completely contained in the 0(5) subgroup of SU(1, 1) x 0(5), 
no explicit representations in terms of 7 had to be constructed to obtain matrix elements 
of the collective variables, relevant for general potential energy surfaces. This suggests 
that the proposed technique can be extended to higher rank algebras, such as the 0(7) 
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orthogonal group, emerging from octupole degrees of freedom in atomic nuclei. 
Now the theoretical framework is set, it is interesting to study to what extend the 
geometrical model can be applied to the collective behaviour of atomic nuclei with 
respect to the recent developments in exotic nuclei. This will be the subject of further 
investigations. 
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Appendix A. Explicit expressions of the 0(5) generators in the 
Cartan-Weyl basis 

The generators Lm and Ou can explicitly be expressed in terms of the collective variables 
and their canonic conjugate momenta according to the definition (jSJ) 

• rr ■ m (A ' 1} 

where (ji^ii2^2|i3^3) are the commonly known Clebsch Gordan coefficients. Taking 

the rotation to the Cartan representation into account ([HD, explicit and relatively simple 

expressions for the generators can be obtained 

% % 
X + = t("27t* i - a-ivr^), Y + = ^(a 2 7T* - ai7r|), 

h h 
— % — % 

X- = — (ai7r* 2 - a^ 2 irT), Y_ = — (a_ 1 7r* 2 - «_ 2 tt* i), 

n a 
— i 

X = ■^■(a 2 vr* 2 + a 1 ir*_ 1 - a_i7r* - a^l), 

— 1 



(A.2) 



— 1 — 1 

— z — z 
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